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1 Introduction 


Recent discovery of the Higgs particle [T] suggests that the framework of the standard 
model (SM) can describe Nature up to the weak scale. On the other hand, we have 
several experimental results which are not explained within it. Representative ones are 
the existence of both small neutrino masses [2] and dark matter (DM) [3] , and also the 
baryon number asymmetry in the Universe |1]. These require some extension of the SM. 

As such an extension, we have an interesting simple model, in which the SM is extended 
by an inert doublet scalar and singlet fermions [5]. It has promising features for the 
simultaneous explanation of the neutrino oscillation data [2] and the observed abundance 
of DM [3] through physics at TeV regions. In fact, if these new helds are assigned odd 
parity of an imposed Z 2 symmetry, the neutrino masses are generated at one-loop level and 
the lightest Z 2 odd particle can behave as DM. The quantitative conditions required for 
their explanation and also other phenomenological aspects have been extensively studied 
in this model and its extended models piTiiHiEiiiniin]. They show that the simultaneous 
explanation of these is possible without causing strong tension with other phenomena like 
lepton flavor violating processes as long as DM is identihed with the lightest neutral 
component of the inert doublet scalar Moreover, in that case, the baryon 

number asymmetry in the Universe could be also successfully explained if rather mild 
mass degeneracy is assumed among the singlet fermions with masses of TeV scales m- 

In this paper, we consider how inflation can be embedded in this framework. CMB 
observations [I31[I1[I5] suggest that the exponential expansion of the Universe should 
occur before the ordinary Big-Bang of the Universe. On the other hand, the analyses of 
them seem to have already ruled out a lot of inflation models proposed by now. Higgs 
inflation is a well-known example which is still alive [16]. This model is characterized 
by the feature such that Higgs potential becomes flat enough for large held regions if the 
Higgs scalar has a large non-minimal coupling with the Ricci scalar. We apply this idea to 
a radiative seesaw model extended by real singlet scalars im. Although the singlet scalars 
are originally introduced with the aim of generating the neutrino masses, it could work as 
inhaton if they are supposed to have a substantial non-minimal coupling with the Ricci 


"‘if the lightest singlet fermion is identified with DM, strong tension appears between the DM abun¬ 
dance and the lepton flavor violating processes [6]. However, it could be resolved by assuming special 
flavor structure [7] or introducing a new interaction [Sj. 
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scalar. In fact, such a coupling of a real singlet scalar has been studied as s-inflation in a 
different context [IB]. Following it, we focus our attention on such a non-minimal coupling 
of the singlet scalars with the Ricci scalar instead of the one of Higgs scalar and others. 
In this case, unitarity problem which appears in the Higgs inflation and other general 
models nnuBU] might be escaped under certain conditions. Moreover, the singlet scalars 
could also play an important role in the generation of the baryon number asymmetry in 
the Universe through non-thermal leptogenesis. We study this issue intensively. 

The following parts of the paper are organized as follows. In section 2, we introduce 
the model studied in this paper and discuss the neutrino mass generation. In section 3, 
a possible inflation scenario in this model is discussed. Production of the baryon number 
asymmetry due to the inflaton decay is studied in detail in section 4. Consistency of the 
DM physics with this scenario is also discussed here. Section 5 is devoted to the summary 
of the paper. 

2 An extension with real singlet scalars 

The radiative seesaw model proposed in [5] is characterized by a scalar quartic coupling 
between the ordinary Higgs doublet cf) and the inert doublet scalar rj. Since 
r] and singlet fermions W are assigned odd parity of the Z 2 symmetry and all the SM 
contents are assigned even parity, the Dirac neutrino mass terms are forbidden at tree 
level. Neutrino masses are generated through a one-loop diagram with W and rj in the 
internal lines. In this mass generation scenario at TeV scales, the above mentioned quartic 
coupling between 0 and rj plays an essential role to explain the small neutrino masses. 

An extension of the model might be done by considering a possibility that this quartic 
coupling is an effective coupling appearing at low energy regions after integrating out 
heavy scalar helds irzi. Such a scenario could be realized by introducing Z 2 odd real 
singlet scalars Sa {a = 1,2)1^ The model is dehned by a part of Lagrangian relevant to 

^One real scalar is enough for the neutrino mass generation and inflation. However, if we consider 
leptogenesis in the model, two real scalars should be introduced at least. We take this minimal version 
here. 
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Fig. 1 The one-loop diagram which contributes neutrino mass generation in the present model. 


the new fields as follows, 


Mk ,,, , Mk 


hakNkV^L + hl,iaVNk + + -Y^kNk 


Q,fc = l *- 

+ m1ct)^(l) -f rn^rj^T] + Ai(0V)^ + ^2{v^vY + + -^4(hV)(0^h) 


E 

a=l,2 


+ ^5'^(0V) + '^Sl{T]^r]) + /ia^ahV + fJ‘*aSa<P^V 


( 1 ) 


where £a is a left-handed doublet lepton. We note that A 5 (? 7 V)^ is allowed under the 
imposed symmetry in general. However, if we assume A 5 = 0 in the original Lagrangian, 
its /3-function is proportional to itself and then A 5 = 0 is stable against the radiative 
correction as long as pa terms are not included in eq. ([I]). On the other hand, if the 
Ha terms are introduced in eq. ([I]) assuming A 5 = 0, the A 5 term appears effectively as 
discussed below. Later, this point will be discussed again in relation to the assignment of 
lepton number to the new fields. 

In this model, neutrino masses are generated through a one-loop diagram which is 
shown in Fig. 1. They are estimated as 




0 

hakhi3kMk{(j))‘^ \ 

2 ^- ^2 - Mk, m,J, 

k=l a 


( 2 ) 


where = mi + (A 3 -b \^){4>)‘^ and (0) = 174 GeV. The function I is defined by 


I{ma,mb,mc) = 


(m^ — m^m^) Inm^ 


+ 


mi In mi 


{ml — miy {m^ — miy {ml — ml){ml — mfy 

ml In ml 1 

{ml - ml){ml - ml)^ {ml - ml){ml - ml)' 


( 3 ) 
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If we suppose ms^ ^ Mn, this formula can be approximated as 




hakhj3k{ft^') 


m% — 

a k=l 


Svr^ 


M^-M^ L 


Mi , Ml ■ 

-i-In —h + 1 

Mi-Mi 


( 4 ) 


„ is idem 

CL ml 


This is equivalent to the neutrino mass formula in the original model if ^ 

Sa 

tihed with the quartic coupling constant As HT]. This correspondence could be directly 
conhrmed in an effective model at energy regions smaller than ms^, which can be derived 
by integrating out Sa- In fact, since the equation of motion for Sa could be approximated 
as Sa — -fhs required terms are derived by using it as 






mt 


. , (5) 

L c TTl Q 

„ ‘^a iJa 

Origin of the smallness of IA 5 I, which is a key to explain the small neutrino masses in 
the original model, is translated to the hierarchy problem between /!„ and ms^ in this 
scenario. We cannot answer the origin of this hierarchy at the present stage and we have 
to leave it for a complete theory at high energy regions. 

For the later study, we show an example of flavor structure of the neutrino Yukawa 
couplings which can explain every neutrino oscillation data in the normal hierarchy case. 
Here, we follow the procedure given in m- For this purpose, we assume that the neutrino 
mass matrix (|4j) takes the following simple form as 



/ 

0 

0 

0 
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/ 

1 

92 

-93 




M’' = 


0 

1 

Qi 


{h\Ai + hlA2)-^ + 


q2 

9 i 

-9293 



( 6 ) 
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where ^ 


‘Si 


IS assume 


‘S 2 


Afc — 


C 0 


and Ak is represented as 




IGeV 

Mfc 


87 r 2 GeV ^ _ 1 
mF 


1 + 


Mi 


Mi 

Mi-Mi ^ Mi, 


{<t>? 


87 r 2 GeV 


At 


(7) 


If we put gi, 2,3 = 1 in eq. ([ 6 ]), the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing 
matrix is found to reduce to the tri-bimaximal form 


UpMNS— 


/ 

2 

P6 

1 

P3 

0 
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0 

0 
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-1 

1 

1 
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^iai 

0 



P6 

P3 
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0 

0 

^ia2 

/ 


( 8 ) 


Ve V3 


®As explained in the later discussion, this assumption is adopted in connection with leptogenesis. 
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Fig. 2 A region in the (92,®) plane allowed by the neutrino oscillation data, which is included inside 
a circle drawn by the black solid line. Other parameters are fixed to satisfy the condition (flOl) . Each 
contour represents 2 (t boundary values of neutrino oscillation parameters |21j . that is, |Am|2| (thick red 
solid and dashed lines), Am^i (thin red solid and dashed lines), sin^ 2023 (green solid and dashed lines), 
sin^ 2012 (blue solid and dashed lines) and sin^ 20i3 (black solid and dashed lines). 


where Majorana phases 01^2 are determined by the phases hi and /ta- If we put ipi = arg(hj) 
and = arg(/ia), they are expressed as 


Ctl — <y53 + a2 — P2 + Pii2-, 


( 9 ) 


where |h.i| is taken as a negligibly small value compared with others, for simplicitylll 
Since one of mass eigenvalues is always zero in this flavor structure, we hnd that the mass 
eigenvalues should satisfy ^ and ~ where 


atm 


t- o ^ I I t fy 

02 02 

and Am^Qi stand for the squared mass differences required by the neutrino oscillation 
analysis for both atmospheric and solar neutrinos mm- 

Since we now know that sin^is takes a non-zero value, we have to consider a flavor 
structure deviated from ^ 1 , 2,3 = 1- For that purpose, we determine values of ^ 1 , 2 , 3 , ^1^2^ 

2 ' ' ’ ’ •s'2 

and h|A 3 ^^^ to realize all the squared mass differences and the mixing angles required 
^2 


^The model is equivalent to the one with two singlet fermions in this case. It should be noted that 
the neutrino oscillation data could be explained as long as only two singlet fermions are introduced. We 
use this setting throughout the following study. 
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Table 1 Examples of the model parameters which satisfy the latter two conditions in eq. (flOl) . If we fix 
a point in a circle of Fig. 2, all the neutrino oscillation data can be explained at 2cr level, consistently. In 
all cases, too, = 10® GeV and = 1.1 are assumed. 

by the neutrino oscillation data through diagonalizing the matrix ([ 6 ]) numerically. This 
analysis shows that the neutrino oscillation parameters can be in the 2 cr range of the 
experimental data if the values of {q 2 , qs) are contained in the region surrounded by a 
circle in Fig. 2. In this figure, the remaining parameters are hxed so as to satisfy 

qi = 0.77, = 6.03 x 10■^^ = 1.02 x 10"^^ (10) 

As long as the model parameters M^, ha ,3 and Ma ,3 are varied by keeping these 
conditions, the neutrino oscillation constraints are automatically fulhlled. In Table 1, 
typical examples obtained by this simple procedure are shown. They include examples 
such that the masses of the singlet fermions are largely different in the case = 
1 TeV, in which they are of O(IO^) GeV and 0(10®) GeV in the cases (a) and (b), 
respectively. 

If Sa does not play any other role than the neutrino mass generation, this extension 
may not be so appealing. However, we can find that the introduction of Sa could add 
favorable features as an inflation model to the radiative seesaw modelj^ Recent Planck 
data suggest that the Higgs inflation scenario could be one of the favored inflation models. 
However, if multi-component scalars like the Higgs doublet scalar are supposed to play a 
role of inflaton, the model could be suffered from the unitarity problem [191 EQ]. Since 
unitarity could be violated in the scattering amplitudes among scalars with non-minimal 

®We have proposed another inflation scenario in the similar context based on somewhat different 
motivation in [22) . The inflaton potential assumed there differs from the present one. As a result, the 
predicted values for the spectral index and the tensor-to-scalar ratio take distinct values from the present 
ones. 
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couplings to the Ricci scalar at a lower scale compared with the inflation scale, new physics 
required for unitarity restoration could jeopardize the flatness of inflaton potential at the 
inflation scale. The situation can be changed in a real singlet inflaton as discussed in [20] , 
In the following part, we consider that only the singlet scalars among scalars in the model 
have non-negligible non-minimal conplings with the Ricci scalar. 


3 Inflation due to the non-minimal coupling 


It has been known that a scalar held coupled with the Ricci scalar can bring about the 
exponential expansion of the Universe [23]. Using this idea to the SM, Higgs inflation 
has been proposed in [16] as a scenario with a realistic inflaton candidate. After that, 
the scenario has been studied from various view points [23]. We apply this idea to the 
singlet scalars in this model but not to the Higgs donblet or the inert donbletE The 
action relevant to the present inflation scenario is given in the Jordan frame as 



a a 


( 11 ) 


where Mpi is the reduced Planck mass and V{Sa) stands for the corresponding part of 
the Sa potential in eq. ([I]). We note that only the singlet scalars are assnmed to have 
non-minimal coupling with the Ricci scalar. 

We take Si as inflaton and other scalars are assnmed to have much smaller values than 
Si during the inflation. In that case, V{Sa) can be approximately expressed as V{Sa) — 
for a sufficiently large value of Si, where the coupling k) of inflaton is abbreviated 
as Ki and HiSf 3> is supposed implicitly. In order to derive the corresponding action 
to eq. flTT]) in the Einstein frame, we use the conformal transformation [T6l 123] 

J = Je, SJ" = 1 + (12) 

As a result of this transformation, we find that it is written as 


Se = d xi/-gE 






a,b=l,2 


K 




(13) 


®The study of Higgs inflation in the inert doublet model can be found in [^. Although the present 
inflation scenario and its prediction are essentially the same as [ISllIH], we note that the inflaton is shown 
to play crucial roles in the neutrino mass generation and the leptogenesis in this model. 










An important feature is the appearance of the mixing between d^Sa and dfj_Sb in the 
second term. We will discuss it later. 

We consider a case in which only one real scalar Si has the non-minimal coupling with 
the Ricci scalar at hrst. In that case, a canonically normalized held y can be introduced 
as 


dx 

dSi 


s? 


i + (ei + 6e?)^ 




1/2 


1 + 


65? 


(14) 


The potential can be expressed by using this y. It is easily seen that the new 


held X coincides with Si at the regions where Si is satished. On the other hand, 


if Si takes a large value such as Si S> Si and y are found from eq. (IT4l) to be related 
as Si oc exp ( —^= -). The potential at this region is found to be almost constant as 






KiMf 


pi 


4(1 + ^ 


4^? 


(15) 


This suggests that y could play a role of slow-rolling inhaton in this region. 
The number of e-foldings induced by the potential can be estimated as 


N = 




dx 




( 1 ) 


R 


( 1 )' 


Q C 2 _ C 2 

^l,end 

4 ’ 


( 16 ) 


where V^'^' = 
are summarized as 


and eq. ffTTl) is used. Slow roll parameters derived from this potential 

[26] 


1 /Ri')'y_4Mp4 



4M2 

?.iiSl 


(17) 


Since the inhation is considered to end at e ~ 1, we have Rf end — v/ which suggests 
that Ri^end could be neglected in eq. flT^ . Thus, the slow roll parameters are found to be 
expressed as £ ~ and 7] ~ —N~^ by using the e-foldings N only. 

The spectrum of density perturbation predicted by the inflation is expressed as 


If we use As = 
Kl ~ 10“®^fiV 


V{k) =As{- 


n. —1 


R 


( 1 ) 


24ir2My 


( 18 ) 


(2.445 ± 0.096) x 10 ® at R = 0.002 Mpc ^ [13], we hnd that the relation 
should be satished at the horizon exit time of the scale fc*. The spectral 
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index Ug and the tensor-to-scalar ratio r are represented by using the slow-roll parameters 



ng = 1 — 6e + 2r], r = IQe. 


(19) 


Using eq. ffT7|) to these formulas, they are found to be determined only by the e-foldings 


N such as Ug ~ 0.968 and r ~ 3.0 x 10 ^ for N = 60. These values coincide with the 


ones estimated from the Planck data well. Although all these results are the same as 
the ones found in the Higgs inflation, the quartic coupling ki is a free parameter in this 
model. It is completely different from the Higgs inflation where the corresponding quartic 
coupling is constrained by the Higgs mass 126 GeV. As a result, we cannot relate weak 
scale physics to the inflation through the observational data of the Universe in this model. 
On the contrary, this fact allows that takes a much smaller value compared with the 
one of the usual Higgs inflation. For example, = O(IO^) can realize both A^ = 60 and 
the observed value of Ag if a very small value such as 0(10“®) is assumed for ki. However, 
as found from the expression of the slow-roll parameters which depend only on N^ the 
predicted values for Ug and r are the same as those of the Higgs inflation. 

Next, we consider the model with two real scalars, which corresponds to the one 
discussed in the previous section. The situation could largely change if multi-scalars have 
couplings with the Ricci scalar. In general, the mixing in eq. flT3|) cannot be resolved by 
any field redehnition and then it is difficult to hnd canonically normalized basis for them. 
An exceptional situation for this could be found for hierarchical couplings such as 3> 1 
and ^ 1^2 1- This condition can be freely imposed on the present model since Si and S 2 

are not related by any symmetry. In that case, the model could behave as a single real 
field model [20] . We can introduce a canonically normalized field y for Si in the same way 
as eq. flTTl) . On the other hand, the S 2 relevant terms in eq. flTdll are strongly suppressed 
as long as ^iSf > is satisfied. In the potential V, the S 2 relevant part can be given 
as 



( 20 ) 


This means that only the y could play a role of slow-roll inflaton also in this case. 

In addition, under the condition ,^ 1,^2 1, the scale of unitarity violation could be 


comparable to the inflation scale [20]. The unitarity violating scattering induced by 


the mixing part in the second term in eq. flTS]) could give the strongest constraint. A 
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simple power counting for the scattering amplitude between Si and S 2 suggests that the 
unitarity violating scale is given by A = However, since the condition \/^i ^2 < \/lT 

is satished, the unitarity violating scale A could be comparable to or larger than the 
inflation scale ^=. Since Si and S 2 are independent helds in the present model, ^2 1 

is possible even if we assume a suitable value of for inflation. Other possible unitarity 
violation induced by other parts such as Ve might be studied through the analysis by 
taking account of the background held dependence. It suggests that the unitarity violation 
scale is comparable to the inhation scale or larger than that. Thus, the hatness of the 
present inhaton potential is reliable throughout the inhation period. Any physics which 
remedies the unitarity violation does not ahect the present inhation scenario. 


4 Non-thermal leptogenesis and dark matter 


4.1 Leptogenesis 

Reheating after inhation is another important problem for this inhation scenario to be 
realistic. If we impose the existence of sufficient thermal relics of the inert doublet DM 
in the present Universe, reheating temperature should be higher than its mass at 
least, which is supposed to be of 0(1) TeV in the present study. Since the allowed decay 
mode for the inhaton Si is limited to Si —?■ the decay width of Si could be 

estimated as T^^ = Applying instantaneous thermalization approximation to this 

process, the reheating temperature could be estimated from the condition ~ T^^ as 

r« ^ 1.6 X 10‘^ (47M77) ^ GeV. (21) 


lOSGeVy 

where 5 ^* = 116 is used. If Tr > is satished and Pr exists in the thermal bath, the 
observed DM abundance could be explained as the relic abundance of this thermal 
The relic abundance is discussed in the next subsection. 

Several leptogenesis scenarios may be considered in this reheating processe depending 
on the lepton number assignment for new ingredients. First, we consider an ordinary 
lepton number assignment to the new helds p and such that L{ri) = 0 and L{Nk) = 1 . 
In this case, if the reheating temperature is higher than a certain bound required for the 


heavy singlet fermion masses 


ij the usual thermal leptogenesis could work. On the other 


discussed in m, this bound for the singlet fermion masses could be relaxed in the radiative seesaw 
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hand, if the reheating temperature is not so high but high enough to thermalize the singlet 
fermions with masses of 0(1) TeV for example, the sufficient baryon number asymmetry 
could be generated through the resonant leptogenesis [28l |29] as long as the masses of the 
singlet fermions are hnely degenerate [m. 

If we hnd that there is another possible assignment of the lepton number such as 
L{ri) = 1 and L{Nk) = 0 |T3 130], we can consider a new leptogenesis scenario allowed 
in this model. It is based on the non-thermal generation of lepton number asymmetry 
through the inflaton decayo Although the quartic coupling is forbidden by this 

lepton number assignment, it could be generated effectively through the lepton number 
violating tri-linear scalar couplings pa at low energy regions as discussed in the previous 
section. Since this coupling violates the lepton number, the decay of inflaton induced 
through this coupling could generate the lepton number asymmetry in the rj sector through 
the interference between tree and one-loop processes. The CP asymmetry expected in 
this decay can be estimated as 


r(^i ^ _ f (Ai ^ 7^10) 


Ih2 


Stt 


In 


rrig^ + ^ 


52. 


+ 


mg, - 


m 


Si 


m 


52 


[m 


Si 


m 


52 




sin2(0i-02), (22) 


where 9a = aig{^a) and la the following study, we assume the maximum 

CP phase | sin2(0i — 6 * 2 )] = 1. 

Since 77 ^ is supposed to be DM, all components of 7] is expected to be lighter than 
Nk. In that case, the generated lepton number asymmetry cannot be transferred from 
the r] sector to the doublet lepton sector through the r] decayo However, it can be 
converted to the lepton sector through 2-2 scattering processes. If this conversion occurs 


model in comparison with the famous Davidson-Ibarra bound in the usual seesaw model m- 

^^The generation of lepton number asymmetry through the inflaton decay has been considered. In m, 

for example, the asymmetry is supposed to be generated by its decay to the heavy right-handed neutrinos 
and their successive decay in the SO(IO) GUT framework. Mass of the decay products is largely different 
from the one in the present scenario. 

^^If a singlet fermion is considered to be DM, rj could decay to the lepton and then the lepton number 
asymmetry in the 77 sector moves to the lepton sector directly through it m- In this case, unfortunately, 
the relic DM abundance could have serious tensions with the lepton flavor violating processes (LFV) [6] . 
However, if we assume a certain flavor structure for neutrino Yukawa couplings, the LFV constraints 
could be satisfied. Detailed analysis for realistic parameters will be presented elsewhere. 
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efficiently without inducing any contradiction with other phenomenological constraints, 
the sphaleron interaction is expected to generate the baryon number asymmetry from this 
lepton number asymmetry. It has already been studied in a different inflation scenario 
[30] . However, since the parameter space of the present model is much simpler than the 
one of that model, we can study the feature of the scenario in a systematic way. We start 
with a brief review for the method of the analysis at first. 

The lepton number asymmetry in the co-moving volume is expressed by using the 
entropy density s as Al^ = jn the doublet lepton sector and as Al^ = ^ 

in the rj sector, respectively. Boltzmann equations which describe the evolution of these 
quantities are given as 


dAYr, 

dz 

dAYL 

dz 


/AW AYi 


sH{Mr 


y 2(7a + 7b) ~ + 2 ( 73 ; + 7^)^^ 


AW 


sH{Mrj 


-2(7a + 7fe) 


/AW AW 


Ur 


i 


n'" J ’ 


(23) 


where we introduce a dimensionless parameter 2 ; which is defined as z = ^. The equilib¬ 
rium values for these are expressed as Y^^{z) = ■:^^z‘^K 2 {z) and Yp = where K 2 {z) 
is the modified Bessel function of the second kind. These equations are derived under the 
assumption such that the existence of S' 1^2 in fhe thermal bath can be neglected. It means 
that both the inverse decay of 81^2 and the scattering containing S' 1^2 in fhe initial and hnal 
states do not contribute to these equations. This requires Tr < ms^ 2 reheating 


temperature Tr. If we take account of eqs. oa and dUD with it, ^ < 


l/^lP ^ 1^2^ 


‘Si 


-T- is found 

’32 


to be imposed. Reaction densities ■ya,b for lepton number conserving scattering processes 
rjr] iai/s and 7 /i pp cause the transition of the lepton number asymmetry between 
the 7 ] sector and the doublet lepton sector. On the other hand, reaction densities ■jx,y for 
lepton number violating scattering processes 7777 —?> 00 and rjp —)■ pcj) control the washout 
of the lepton number asymmetry in the rj sector. Formulas of these reaction densities are 
summarized in Appendix A. If we note that 77 , 5 j(Tr) = and P 5 ^(Tr) = are 


satisfied for the number density and the energy density of Si at Tr under the assumption 


of instantaneous thermalization, A10 (Tr) = 


could be obtained by using eq. 

We use it and AYl{TP) = 0 as the initial values for this analysis. The lepton number 


asymmetry AYl obtained at the weak scale as the solution of these Boltzmann equations 
is converted to the baryon number asymmetry through the sphaleron processes. In the 
present model, the resulting baryon number asymmetry could be estimated by using the 
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Fig. 3 Relevant reaction rates and solutions of the Boltzmann equations for the cases (a) (upper panels) 
and (b) (lower panels). They correspond to the cases with single fermion masses of O(IO^) GeV and 
0(10®) GeV. In the left-hand panels, the ratio of reaction rate T to the Hubble parameter H for the 
relevant process is plotted as functions of z. In the right-hand panels, the evolution of AV, and AY^ are 
plotted as functions of z. The lepton number asymmetry required to explain the observational results is 
shown by the horizontal black dashed line. 


solution of eq. fl2^ 



Yb = - 


7 

19 


^Yl{zew)- 


(24) 


Now we show the results of the numerical analysis of the baryon number asymmetry 
generated through the scenario described above. Since a factor is included 

in the reduced cross section da^b given in Appendix A, we have to determine the flavor 
structure of neutrino Yukawa couplings for the realistic analysis. We adopt the results 
obtained in Fig. 2 for it and then they are hxed as 


{hh))22 — (1 + (hh^)33 — (1 + Q'2 + (25) 

find a relation B = —^{B — L) at the weak scale from the chemical equilibrium condition [T7]. 
Since B — L is conserved under the sphaleron interaction, we note that eq. (1^ should be regarded as 
the Boltzmann equations for it but not for L. 
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where we could choose the values of ^ 2,3 in the allowed region shown in Fig. 2 for qi = 0.77. 
In the following analysis, we use q 2 = —0.4 and ^3 = 2.1 which fix the neutrino oscillation 
parameters as Am |2 = 2.51 • 10“^ eV^, Amh = 7.63 ■ 10“^ eV^, sin^ 2023 = 0.977, 
sin^ 2012 = 0.863 and sin^ 20 i 3 = 0.097. In Fig. 3, we show the solutions obtained for the 
model parameters given as the cases (a) and (b) in Table 1. For each case, the ratio of 
the relevant reaction rate to the Hubble parameter is plotted as a function of 2 ; in the left 
panels. The evolution of AY^j and is shown in the right panels, where AYl{zew) 
required for the observed baryon number asymmetry is implicated by the horizontal dotted 
lines. We hnd that the sufficient amount of baryon number asymmetry is generated for 
parameters which are consistent with the neutrino oscillation data. The obtained values 
are listed in the last column of Table 1. 

In Fig. 3, we find completely different behavior in the transition of the lepton number 
asymmetry between the two cases. The difference comes from whether the lepton number 
conserving scatterings could be in the thermal equilibrium or not. It is determined de¬ 
pending on both values of the neutrino Yukawa couplings and singlet fermion masses. In 

this model, they are constrained by the neutrino oscillation data together with the value 

2 _ 

of the effective coupling A 5 (~ -^) as found in eq. (ITOl) . If the lepton number conserving 

^S2 

scatterings could be expected in the thermal equilibrium for Tr ::§> M^, the situation 
AYl — AYrj is realized during that period. It is found in the case (a). The final value 
of AYl at the electroweak scale is determined depending on the time when they leave 
the equilibrium. On the other hand, if we suppose that these processes could never be 
in the thermal equilibrium, we find that Tr <C Mk should be satisfied and the relation 
AYl — AYr, cannot be kept during any substantial period. This corresponds to the case 
(b). The final value of AYl in this case is hxed mainly by the strength of the lepton 
number conserving scattering processes at Tr. 

It is useful to clarify the parameter dependence of the generated baryon number asym¬ 
metry Yr in order to understand this scenario. In the case (a), Yr is stable against the 
change of 0 . 2,3 and M 2,3 which satisfies the conditions given in eq. flTOl) as long as A 5 is 
kept to be a constant value. This feature is considered to be brought about since the 
decoupling temperature of the lepton number conserving scatterings is not affected by 
this variation substantially. On the other hand, Yr is deeply dependent on the values of 
Pi and p 2 - The former one determines the initial value of AY,^ and the latter one controls 


15 


10 '® 
10 '® 
10 "'“ 




10 *' 


10 *' 


10 *' 




+ 

X 

+ 

X 

+ 


A 


A 

A 

X 

X 





A 


|?.5l=1.0x10*'' 

+ 





■ I1.5h1.5x10*'' 

X 





hjhZ.OxlO*'' 

A 






10*® lO*'' 

I^l1l/ms, 



Fig. 4 Left: The dependence of IYbI on \/[A7f) which are relevant to the initial value of the 

lepton number asymmetry and the washout of the lepton number asymmetry, respectively, msj = 
is assumed. Right: The dependence of |Yb| on A = — 1 for several values of [As]. = 2 • 10“® 

is assumed. In both panels, other parameters are fixed to satisfy the neutrino oscillation data assuming 
msi = 10® GeV and = 1 TeV as discussed in the part relevant to eq. (ITOll . 


the washout of Examples of this dependence could be seen are through the left 

panel of Fig. 4. Smaller |/ii| makes the reheating temperature Tr lower and then reduces 
the initial lepton number asymmetry AYj^(Tr). Larger |pi| makes Tr higher and Tr could 
take a near value to fhe present setting. This makes the washout of at the 

neighborhood of Tr be enhanced due to the tail effect of the s-channel resonance. It can 
be seen at a small 2 ; region in the figures of These could explain the reason why \Yb\ 
is a convex function of Since larger IA 5 I makes the washout effect larger as found 
from the formulas of ax,y given in Appendix A, IA 5 I is expected to have an upper bound. 
The left panel of Fig. 4 shows that IA 5 I should be smaller than 2 ■ 10“® to guarantee the 
sufficient amount of IYbI hi this case. Degeneracy between ms^ and ms 2 is also crucial 
to fix the value of the CP asymmetry and then the initial lepton number asymmetry 
AYy{TR). In the right panel of Fig. 4, \Yb\ is plotted as a function of A(= — 1). If 

the mass difference becomes smaller, the second term in eq. fl 2 ^ which comes from the 
self-energy diagram is enhanced and then the initial value of AY^{Tr) becomes larger. 
However, as shown in these examples, the fine degeneracy is not required as long as A 5 
takes a smaller value as mentioned above. It is a distinctive point from the ordinary TeV 
scale thermal leptogenesis. This comes from the non-thermal origin due to the inflation 
decay which could prepare a sufficient amount of asymmetry. 
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In the case (b), the situation is completely different from the case (a). Yb is largely 
dependent on the setting of h 2,3 and M 2^3 for a fixed value of A 5 . Unless is much 
larger than Tr, the rate of the lepton number conserving scatterings could be enhanced 
to be almost in the thermal equilibrium at the neighborhood of Tr. In that case, since 
the washout processes are in the equilibrium due to an assumed large value of IA 5 I, 
decreases steeply and the transferred AYl follows AY^j not to reach a substantial value. 
To escape this situation, it seems to be necessary for to be larger than Tr by an order 
of magnitude at least. These features clarify how both and play crucial roles in 
this scenario. Anyway, these studies suggest that the present non-thermal leptogenesis 
could give us a successful scenario for the generation of the baryon number asymmetry. 
It is completely consistent with the neutrino mass generation which could explain all the 
neutrino oscillation data. 


4.2 Dark matter 


In this subsection, we discuss the connection between this leptogenesis scenario and DM 
physics under the assumption that the lightest neutral scalar r]\ is DM. At first, we 
discuss the relic abundance of r]\. In the mass range of discussed in this paper, its 
abundance is known to be determined by the couplings A 3^4 in eq. o jTTiiia, which are 
completely irrelevant to the analysis of other phenomena studied hereo As discussed in 
ED], the lepton number asymmetry kept in the 77 sector cannot play any role in the DM 
abundance. The asymmetry in this sector disappears through the effective coupling A 5 
after the electroweak symmetry breaking. Thus, the required relic abundance should be 
realized as thermal relics for suitable values of A 34 . The estimation of its relic abundance 
for Mri = 600 GeV, 1 and 3 TeV is shown in Fig. 5, which is obtained by using the ordinary 
method whose detail can be found in [30]. Since the required abundance is displayed by 
the horizontal dotted line in this hgure, we hnd that appropriate values of A 3^4 could make 
the thermal t^r be a favorable DM candidate naturally. 

Next, we discuss the consistency between this leptogenesis scenario and the r^R DM. 

2 

Since the effective coupling A 5 (~ -^) takes a small value, the real and imaginary parts 

S2 

of the neutral component of t] have almost degenerate masses whose difference is 


note that the required relic abundance cannot be explained for < 530 GeV where only the 
gauge interaction reduces it below the required value [12) . 
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Fig. 5 Thermal relic abundance for = 600 GeV, 1 TeV and 3 TeV from the left to the right, 
respectively. It is estimated by taking account of the coannihilation between the components of 7] which 
is controlled by the couplings Aa^. 


expressed a 



6 = Af„o — ibf„o 

'll 'In 




On the other hand, rj^ and rj^ have a weak gauge interaction such as 


( 26 ) 


^int = ^Z^{rindn,r]i - r]idn,T]n). 


( 27 ) 


This interaction induces the spin-independent inelastic scattering rjuN rjiN with a 
nucleus N mediated by a exchange in which the nucleus is not excited as long as the 
mass difference 6 is sufficiently small. This reaction can contribute to the direct search of 
DM. The r/^-nucleon cross section for this inelastic scattering can be estimated as 

c^iinei = ^ 7.44 X 10 - 3 ^ cm^ ( 28 ) 

27r/i^ 


where fin is the reduced mass of this r^^-nucleon system and 6 is assumed. If we 

apply the bound from recent DM direct search experiments to this cross section, we could 
approximately derive a useful constraint on this leptogenesis scenario. 

As briefly described the in Appendix B, the present bound on cr° estimated for the 
elastic scattering (5 = 0 ) might be translated to the one for the inelastic scattering (5 7 ^ 0 ) 
as 


{v-VeT (■v + iier 

dv I e ^ — e ^0 


A Q ^'^min(<^ —0 ) 

^ _ 

n,inel n,el 


(v-ve)'^ _(v±ve)^ 

dv I e ^ — e ^0 




( 29 ) 


77 ^ and 77 ° could be mass eigenstates for a real A 5 . Since we assume that 7 x 2 is real and 
is satisfied, the effective coupling A 5 can be treated as a real parameter. 


Iml" ^ IM2I" 

mi; rric 
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Fig. 6 Constraint on the effective coupling A 5 , which is derived from the DM direct detection experiments. 
Each line represents (t° for A bound for < 7 ° at mDM= 3 TeV, 1 TeV and 600 GeV given 

by LUX is used. 

where Umin is the minimum velocity required to induce this scattering. This cr° 
obtained by using the present bound on cr° should be larger than the value given in 
eq. fl^ since DM has not been detected in any direct detection experiments [32l[33]. We 
also note that 77 ^ cannot be detected in any direct search experiments unless Umin < Wesc 
is satisfied for the local escape velocity Uesc from our Galaxy. It is estimated as 498 
km/s< Uesc <608 km/s [3l] and its medium value 544 km/s is used as Uesc in fhis analysis. 
In Fig. 6 , we plot corresponding to it as a function of IA 5 I for typical values of the 

recoil energy of nucleus Eji by using the bound of cr° given by LUX [33] as in eq. (12^ 
Since the endpoint of each line represents the occurrence of Umin = Vesc, the scattering is 
kinematically forbidden at the larger IA 5 I region than it and then such a region of IA 5 I is 
allowed from the result of the present direct DM search. Unless the escape velocity takes 
much larger values than the one used here, this analysis unfortunately suggests that this 
DM with the mass < 6 TeV is difficult to be detected through this inelastic scattering 
process in the direct search experiments. 

If we take account of the relation IA 5 I ~ in this model, we might roughly read off 

Since the differential WIMP-nucleus cross section in the case of scattering leading to a heavier WIMP 
depends on the recoil energy very differently from that associated with the standard elastic scattering 
and since in the present case the nucleon cross section depends on the WIMP velocity, clearly we cannot 
derive an exact bound by using the experimental result for the elastic scattering. 


19 

















the condition required from this figure as follows, 

{ (1 — 2) X 10“® for = 3 TeV, 

(5 - 7) X 10-» for M,o = 1 TeV, (30) 

(3 - 4) X 10-® for = 600 GeV, 

although it has non-negligible dependence on Er and Vesc- We should note that IA 5 I 
cannot take small values freely even if the neutrino mass constraints have only to be 
satisfied. In the case (b), this constraint is found to be consistent with the parameters 
used in the above analysis of the baryon number asymmetry. On the other hand, in 
the case (a), one might consider from Fig. 6 that this constraint is not satisfied and the 
scenario is inconsistent. However, it may be appropriate to judge that the consistency of 
the parameters used here is marginal if we take account of several uncertainties included 
in the estimation of the bound on IA 5 I. We also note that the situation could be improved 
if we suppose hne mass degeneracy between Si and S 2 such as A = 10“® and a favorable 
value for such as 2 • 10“^. As we find in the right-panel of Fig. 5, these could enhance 
the initial asymmetry of the lepton number produced through the inflaton decay. In this 
case, for example, IFbI = 1.2 ■ 10“^° can be obtained for IA 5 I = 3.5 ■ 10“® at = 1 TeV. 

5 Summary 

We have proposed an extension of the radiative neutrino mass model with Z 2 odd real 
singlet scalars, which give a seed of lepton number violation to allow Majorana neutrino 
mass generation at one-loop level. If they have hierarchical non-minimal couplings with 
the Ricci scalar such as iaS^R with ::§> 1 ::§> ^ 2 , the Si potential at large held regions 
is so hat that sufhcient inhation could be induced. Although both the scalar spectral index 
and the tensor-to-scalar ratio take favorable values just as the Higgs inhation, the model 
can evade from the unitarity problem diherently from the ordinary Higgs inhation. Since 
the unitarity violating scale could be similar to or larger than the inhation scale, the 
hatness of the inhaton potential is never disturbed by new physics which restores the 
unitarity of the model. 

The decay of inhaton could non-thermally produce the lepton number asymmetry 
in the rj sector through the reheating processes. Although its decay cannot yield this 
asymmetry directly in the lepton sector, the lepton number conserving scatterings could 
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convert a sufficient amount of asymmetry from the rj sector to the lepton sector. This 
lepton number asymmetry is transferred to the baryon number asymmetry through the 
sphaleron interaction. Parameters relevant to this leptogenesis are constrained by the 
neutrino oscillation data and the DM direct search experiments. We have shown that the 
sufficient baryon number could be generated consistently with these constraints under 
suitable conditions. Moreover, the DM abundance could be hxed by the parameters 
which are irrelevant to all of these as long as is supposed to be DM. Although the 
model considered here is very simple, it can compactly explain problems in the SM such 
as the neutrino mass generation, the DM origin and its abundance, the inflation and the 
baryon number asymmetry in the Universe. They are closely related each other through 
the radiative mechanism for the neutrino mass generation. 
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Appendix A 


We summarize the formulas for the reaction densities which contribute to the Boltzmann 
equations used in this analysis. We introduce dimensionless variables 


X = 


A/2’ 


— 


Ml 

M2’ 


ba = 


m 


Sa 


Mf 


Cn. - 


Ih'Q 


M 2 ’ 


(31) 


where s is the squared center of mass energy. The reaction density for the scattering 
processes is expressed as 


7(06 -)> ij) = 


T 


647r^ 


ds a{s)y/sKi 


T 


(32) 


where (j(s) is the reduced cross section and Ki{z) is the modihed Bessel function of the 
second kind. The lower bound of integration is dehned as Smin = max[(ma + mM, ij^i + 
m,) 2 ]. 

The reduced cross section &a,b for the lepton number conserving scattering processes 
are given as 

^ '' ak{x‘^ — 4a;)^/2 


da (a;) = 


27r 




k=l 


akX + (ofc - 1)2 
a; + (a;2 — 4a;) ^^2 _|_ 2ak — 2 

X + 2ak — 2 \^a; — (a;2 — 4a;)^/2 _|_ 2ak — 2 

Re[(/i/it)2.]y^aiaJ j 2x + 3ai + — 4 ^ f x + (a;2 — 4 a ;)^/2 _|_ 2a 


In 


X Qji Qjj 


E 

i>j 

2x + ai + 3aj — 4 , /x + (x2 — AxV^'^ + 2a, 

In ' 


Qj — tti 


Qj'i Qj'j 


X — (x2 — 4x)^/2 _|_ 2aj — 2 
-2 


X — (x2 — 4x)^/2 _|_ 2aj 


(33) 


for rjr] iai /3 and 
db(x) = 


1 (x — 1)^ 
27r x2 


) 2 
kk 


X 


(x — 1)2 
^Re[(Mt) 


1 + 


,k=l 

2(x + Ofc - 2) 


X 


+ 


2x 


(x — 1)" 

X ixttk — l Dk{x) ' 2Dk{xy 

x(x + Ofc - 2)' 


Dkix) 


In 


XQjfi 


2 1 

ijl 


i>j 


X X (x — 1)2 

+ + 


X {Di{x) Dj{x) Di{x)Dj{x) 


+ 


X 


(x — 1)^ 
x 2 

(x — 1)^ 


2(x + ai —2) x + Oj —2\ x(x + a, — 2) 

In 


ttj — tti 


Dj{x) 


Xttj — 1 


2(x + ttj — 2) X + Oj — 2^ ^ x(x + ttj — 2) 


a,' — ttn 


DAx) 


Xttj — 1 


(34) 
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for III these formulas, we use the following dehnition for convenience: 

2 


X CL}- 


Dk{x) (x - afc)2 + 0^4’ 


dk = 


647r^ 


5^ \Kk\^ 


\a=e,fi,T 


1 - — 


( 35 ) 


If we take account of the assumption the reduced cross section 4^ of 

777. o 777( o ’ ^ 

bl b 2 

the lepton number violating scattering processes could be approximately represented as 


4(x) ~ 


1 


2 1 P 2 +1 

+ — In 


TT (x 3 (x - 4))V2 \P^-l ' P2-IJ ' 


ay[x) ~ ^ 


2 1 fx - 1)^ 1 1 X - 62 , Q2 + 1 


+ 


+ -- 


In ■ 


(36) 


TT [(X - 1)2 Q2 _ X 2x2 D2{x) X D^[x) ^2 “ 1. 

for rjrj (fycf) and ?70l —?■ rf'cj), respectively. In these formulas we use the dehnition such as 


1 1 1 / \ 2 / 1 \ 2 
1 1 7 


Da{x) [x-haY + hlda 


da = 


647r2 \ba 


1 - 


D _ 2(1-4) -X ^ _ , , 2(1-x4) 

— r/ ,i\n/o5 — -L“r 


[x(x-4)]V2 ’ 


(x — 1)2 


(37) 


Appendix B 

We consider the direct DM detection through the inelastic scattering with the target 
composed of the nucleus with the atomic number Z and the mass number A [35l [36] . Its 
differential detection rate per unit target mass is 


dR 


= W 


Pdm 


d^v vf{v,Ve 


da 


(38) 


dEn m,DM J ' dEji 

where Nt is a number of target nuclei per unit mass and Er is the recoil energy of nucleus. 
The DM velocity distribution in the rest frame of detector may be taken as a Maxwell- 
Boltzmann distribution f{v,Ve) = exp with Vo = 220 km/s. We take 

into account the motion of the Sun and the Earth by using whose magnitude changes 


1.05+ 0.07 cos 


as Ve = Vq 
can occur is estimated as 


1 yr 


^min 


The minimum velocity for which the scattering 


\ rrir J 


(39) 


y/2m]\[Eii 

where uin is the mass of a target nucleus and nir is the reduced mass of DM-nucleus 
system nir = . The differential cross section for spin independent interaction 

•7 ' m-DM+TTiAT dh,T, 
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is expressed by using the DM-nucleon cross section cr° at zero momentum transfer as 


da _ rriN [Zfp + {A - Z)fnf 




dEn /2 

where E{Eji) is a form factor of the nucleus. If we substitute this in eq. 
ential rate can be represented by using eq. fHOj) as 


where is written as 


^ [ d^V - f{v,Ve), 


mjvPDM |Z/p +[A- Z)f„f 


D„ = F\EnY 

2p;mDM /„ 


(40) 
the differ- 

(41) 


We note that Hat takes a hxed value as long as the same target is used. 

The remaining part depends on the sub-process and it can be expressed as 


a'^ [ dv ^ f{v,Ve) =a^ ^ 




y/nVoVe 


_(u —Ue)_ {v + ver 

dv I e ^ — e ^0 


(42) 


(43) 


This depends on whether the scattering occurs elastically (5 = 0) or inelastically (5 7 ^ 0). 
If we interpret the present direct detection results based on these scattering processes, the 
present bound on the elastic scattering cross section cr° gj can be translated to the bound 
on the inelastic scattering cross section jj^gj through 


r'^esc ( (u —ue)^ _ + \ P'^esc / (v — ve)^ _ (‘u-\-ve)^ 

f^n.inel / i 6 ”0 - 6 ”o = gj / dv i 6 "o - e "o 

y J d'V^in{i=0) \ 

(44) 

We use this relation to constrain the allowed values of 6 , which makes us possible to hnd 
the lower bound of the effective coupling [Asj. 
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